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Abstract 

The phase diagram for staggered fermions is discussed in the context of the staggered chiral 
Lagrangian, extending previous work on the subject. When the discretization errors are significant, 
there may be an Aoki-like phase for staggered fermions, where the remnant SO{4) taste symmetry 
is broken down to 50(3). We solve explicitly for the mass spectrum in the 3-flavor degenerate 
mass case and discuss qualitatively the 2 + 1-flavor case. From numerical results we find that 
current simulations are outside the staggered Aoki phase. As for near-future simulations with 
more improved versions of the staggered action, it seems unlikely that these will be in the Aoki 
phase for any realistic value of the quark mass, although the evidence is not conclusive. 

PACS numbers: 11.15.Ha,11.30.Qc,12.39.Fe 
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I. INTRODUCTION 



Current simulations with dynamical staggered quarks have the benefits of both good 
chiral properties at finite lattice spacing and being computationally inexpensive [1]. The 
price for these features is that the quark doubling is not completely removed, and there are 
four "tastes" of quarks for each staggered flavor. Additionally, there are large 0(a 2 ) scaling 
violations that arise and must be taken into account. Improved forms of the staggered 
action (e.g., "Asqtad" staggered quarks) reduce these scaling violations but they are still 
not negligible 

These scaling violations cause mass differences among the various tastes of light pseu- 



doscalar mesons that result from the four quark tastes. In principle, these break the contin- 

n 

uum SU(4) taste symmetry down to the lattice subgroup |7[, however, at 0(a 2 ) there is a 
remnant £0(4) taste symmetry which remains [8]. For each meson flavor, there are 15 tastes 
which fall into representations of SO (A): pseudoscalar (1), axial (4), tensor (6) and vector 
(4); there is also the £77(4) singlet for a total of 16. This approximate £0(4) symmetry is 
broken down to the lattice subgroup at 0(a A ). These splittings may give rise to a negative 
mass-squared of a meson, since there is no reason to assume the splittings are positive. This 
signals an instability of the vacuum, and thus a meson condensate may form. 

The idea that there could be this sort of spontaneous symmetry breaking in the context 
of the flavor symmetry of Wilson fermions was first proposed by Aoki P|. Lee and Sharpe 
showed that this could also occur for one flavor of staggered fermions as well ^. In this case, 
they assumed one of these splittings to be negative and so the vector-taste pion condensed. 
All current formulations of staggered fermions find that these mass splittings are all positive, 
and thus this is not a likely scenario. In the case of multiple flavors, additional parameters 
which mix the axial and vector taste flavor-neutrals are negative in current simulations, so 



ll| than originally thought. The 



a possible "staggered- Aoki" phase may be more likely jlO . 
phase structure of lattice theories has been studied for a variety of other lattice formulations 
recently Qll3, 3 3- 

As we approach the broken phase, the quark mass m becomes comparable to the dis- 
cretization errors; in lattice units, am ~ a 3 AQ CD . It is then that ^i meson may be negative 
and the £0(4) taste symmetry is broken. The breaking is £0(4) down to £0(3), giving 
rise to three approximate Goldstone bosons. We expect a problem to arise in simulations if 
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either very small quark masses or large lattice spacings are used. It seems as though current 
simulations are not in this region of the phase diagram, as we will discuss. Although it is 
not completely ruled out, it appears as though further improvements to the staggered action 
will also be safely away from the staggered- Aoki phase. 

As pointed out in Ref. ^|, we are not breaking an exact symmetry. The SO (4) symmetry 
will be broken down to the lattice subgroup at higher order, so in this way it is unlike 
the Aoki phase for Wilson fermions, where the flavor symmetry is exact. Thus, the three 
Goldstone bosons actually have masses of 0(a 4 ). Terms which arise at 0(a 4 ) are subleading, 
however, and will only give rise to small corrections to our results here. For small enough 
lattice spacing they cannot change drastically the phase structure here. 

The goal of this paper is to study the features of the particle spectrum in the broken 
phase as compared to the unbroken phase. We will do this first for the case where all three 
quark masses are degenerate, the 3-flavor case, which is simple and exactly solvable. We 
will then look empirically at the 2 + 1-flavor case, where the up and down quark masses are 
degenerate but different than the strange quark mass. This case is significantly more difficult 
than the 3-flavor case — we cannot analytically solve for the condensate in this case — and we 
will simply make some general comments on the critical point and how it relates to values of 
the parameters in current and future simulations. We will find that current simulations are 
not in the broken phase, and more improved versions of staggered fermions are likely not 
going to be in the broken phase, although we note that we cannot make exact statements 
about the more improved versions without performing simulations. 

The approach we follow in this paper is similar to that of Refs. 0, 12 , Q, 13 13 > usm g 
the chiral Lagrangian. This approach allows us to determine the mass spectrum of all the 
light pseudoscalar mesons rather simply. For the region we are discussing, m ~ a 2 , the 
chiral expansion is quite natural, and we can use this to determine the meson masses with 
the dependence on all parameters explicit. We are not performing a complete analysis of 
the phase diagram, however. We will see that the phase diagram will turn out to be more 
complex than studied in this paper, but, although this will change certain features such as 
the mass spectrum, it will not change some central results of the analysis, as we will see. 

This paper is structured as follows. In Sec.|Tl]we layout the necessary tools for the analysis: 
the staggered chiral Lagrangian including the 0(a 2 ) taste-violating potential. Next, we 
determine the mass spectrum for the case of 3 degenerate flavors of quarks in Sec. IIHI In 
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Sec. IIVI we discuss the 2 + 1-flavor case. Here we will not solve for the mass spectrum, 
but make some qualitative statements regarding this case. We will examine the features 
of the critical point with current simulations as a function of the quark masses in Sec. 
Additionally we will look at how the point shifts with more improved staggered fermion 
formulations, and will see that indeed the problem becomes less severe in these cases. We 
finish up with some conclusions in Sec. IVII 



111- 



II. THE STAGGERED CHIRAL LAGRANGIAN 

The starting point of our analysis is the SXPT Lagrangian for 3 flavors of quarks jlO . 
The Lagrangian is written in terms of the field E = exp(i$//), a 12 x 12 matrix, with $ 
given by: 

/ TT ^+ TS + \ 



V 



u 


7T + 


K+ 




D 


K° 


K~ 


K° 


S 



(1) 



/ 



where the elements shown are each 4x4 matrices, linear combinations of the Hermitian 
generators 

T a = {£5, Zn, £/}■ (2) 

In other words, U = J2 a U a T a , K° = J2 a Ka^a, etc. In Euclidean space, the gamma matrices 
£ M are Hermitian, and we use the notations = [fi < v in Eq. (J2J], £ M 5 = and 
£1 = / is the 4x4 identity matrix. Under the chiral SU (12) l x SU(12)r symmetry, 
E — > LY,B) . The components of the diagonal (flavor-neutral) elements (U a , D a , S a ) are real; 
while the other (flavor-charged) fields are complex (71"+, K®, etc.), such that $ is Hermitian. 
The Lagrangian is given by 

C = J - Tr(a M Ea M Et) - -af Tr(A^E + ME*) + —±(Uj + Dj + Sj) 2 + a 2 V, (3) 
o 4 o 

where \i is a constant with dimensions of mass, / is the tree-level pion decay constant, the 
/71q term includes the 3 flavor-neutral taste-singlet fields (m — ► 00 at the end to decouple 
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the taste-singlet rfj), and V = U + W is the taste-symmetry breaking potential given by 

-U = CiTr^f^f^) 

+ C '3^[ Tr ^ 3)S ^ 3)S )+ /i - C -] 
v 

V 

+Ce ^X^S^), ( 4 ) 
-W = C 2y ^[Tr(e ) S)Tr(e ) S) + / i .c] 

V 

+C2a\ E^^^) Tr (^ )S ) + ^1 
v 

+^^[Tr(^ 3) S)Tr(^St)] 

V 

+C5A^[Tr(ei 3 5 ) S)Tr(^ ) St)] . (5) 

(3) 

The in V are block-diagonal 12 x 12 matrices 



t(3) 
SB 



^ B o o 

£ B o 

V o o y 



with B G {5, /i5, yUZ/(/i < u),fi, I}. 

The mass matrix, A4, is the 12 x 12 matrix 

/ 



M 



\ 



\ 



J 



(6) 



(7) 



m u I 
m d I 
m s I 

We will in later sections take either all three quark masses to be degenerate or m u = rrid ^ 



As noted in Refs. 
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11], this potential, although breaking the taste symmetry at 0(a 2 ), 



has an accidental 50(4) symmetry. This implies a degeneracy in the masses among different 
tastes of a given flavor meson, which is seen in the tree-level masses of the pseudoscalar 
mesons. We can classify these mesons into irreducible representations of SO (4). The mass 
for the meson M (composed of quarks a and b) with taste B, is given at tree-level by 



rn 



m b = A* ( m a + rn b ) + a 2 A B , 



(8) 



with 





= A P 


= 


A(£„0 


= A, 


= _ (C + 3C, + Ca + 3C fi 




= A T 


= ^(2C 3 + 2C 4 + 4C 6 ) 


A(U 


= Ay 


= ^ (Ci + C 3 + 3C 4 + 3C 6 


A(C/) 


= Aj 


= ^(4C 3 + 4C 4 ) . 



(9) 

As in Ref. we do not include the ml terms or the terms from U ' in these masses. Those 
terms only affect the flavor-neutral mesons and allow for mixings in the basis of Eq. the 
"flavor basis." We can expand IA ' to quadratic order to find these mixings take the same 
form as the singlet term in the Lagrangian: 



taste — vector 
taste — axial , 



(10) 
(11) 



where we have defined the parameters 



V(A) 



16 \r 



c. 



5V(A) 



] ■ 



(12) 



This requires us to resum the flavor-neutral propagators in the three cases: taste-vector, 
taste-axial and taste-singlet [h]. Performing this resummation is equivalent to writing 
everything in terms of the "physical basis" of our mesons, the 7r°, t] and rf for each of 
the three tastes here, which are the true eigenstates of the mass matrix. These are the 
additional parameters discussed in the introduction which may cause the vector or axial rj' 
mass-squared to be negative, thereby introducing a staggered Aoki phase similar to that 
discussed in Ref. gj]. 

The mass eigenstates of the flavor-neutral mass matrix, for non-degenerate flavors, are 
quite complicated, but the completely general case is not needed. The two cases relevant to 
current simulations are the fully degenerate case (m u = = m s = m), which we refer to 
as the "3-flavor" case, and the case where two flavors are degenerate (m u = rrid = m ^ m s ), 
or the "2 + 1" case. We will examine the 3-flavor case in Sec. IHII and discuss the 2 + 1-flavor 
case in Sec. IIVI 



III. MASS SPECTRUM IN THE BROKEN PHASE FOR 3 DEGENERATE FLA- 
VORS 



We first discuss the simplest case, where all three quarks are degenerate in mass, 
m u = ma — m s = m - First we need the eigenvalues of the full flavor-neutral mass ma- 
trix. The vector, axial and singlet tastes are the only tastes which mix and therefore require 
rediagonalizing the mass matrix. The singlet-taste masses do not involve the new param- 
eters given in Eq. ([12)1 . and so we will not focus on their masses. Also, for the singlet, we 
take tuq —>■ oo on theoretical grounds, but there is no similar reaso ning to dictate the values 
of 5y or 5' A ; they must be determined from fits to lattice data jl^l. 

For the 3-flavor case, the mass eigenvalues for the vector (axial) tastes are given by 



rn 



2 2 

V(A) r" 
2 _ o 2 \ , 3 „2,;' 



We have taken into account the transition from four to one tastes per flavor with the addi- 
tional factor of 1/4 in the a 2 S' v ^ term both in Eq. (fTHl) and below. More discussion on how 
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to account for this transition can be found in Refs. [1] 

Lee and Sharpe's original discussion in the one-flavor case focused on the possibility that 
one of the splittings in Eq. may be negative and cause a meson mass-squared to go 
negative for certain values of the quark masses. This signals a vacuum instability and the 
formation of a meson condensate. As noted in Ref. Ja|, all current simulations with staggered 
quarks show that Ab > for all 5^,0]. On the other hand, as is clear from Eq. ()13j) . 
there is another possibility for a meson mass-squared to become negative, when one of the 
"hairpin" parameters, S' v , A y is negative. 1 In current fits to MILC simulations, both S' v and 



5' A are negative 16|, ^ | , so we can see it is more probable to encounter an Aoki-like phase 
with staggered fermions in this way. Current simulations are outside the broken phase p[ H| . 
However, there is the question of how close current simulations are to this broken phase and 
how improving the action more changes this. 



1 The term "hairpin" comes from the fact that these two-point interactions are disconnected at the quark 
level. See Refs. for details. 
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It is straightforward to see that the vector (axial) rf mass-squared is negative if 

4 

a2 $V(A) < a %(A),crit = ~3 [ 2 ^ m + a2A V(A)] ■ (14) 

Aa < Ay in simulations, so Eq. (|T3|) implies that |5^ crit | < | ^v'crit I ■ Thus, 5' A is more likely 
to be less than S' AcTit than for the vector taste. We will assume for the following analysis 
that this is the case, although the analysis for the vector mesons would be analogous. 
To find the vacuum state of the theory, we must minimize the potential 

W = --fj,mf 2 Tr{E + Tt) + a 2 U + a 2 W , (15) 

where we have substituted M. = diag(m/, ml, ml). 

This calculation is most simply done in the physical basis, where everything is written in 
terms of n°, rj and rf instead of the flavor-basis mesons U, D and S. In the 3-flavor case, 
these two bases are related by 

A = 4= (Ub - D B ) , 

Vb = -^=(U b + D b -2Sb) , (16) 
V6 

rf B = -±=(U B + D B + S B ) , 

with B any of the 16 tastes. In the 2 + 1-flavor case, or any case where some masses are 
non-degenerate, this relationship does not hold for all tastes. 2 For the taste singlets, Eq. (|16|) 
is valid in the mo — > oo limit for any quark masses. For the tensor and pseudoscalar tastes, 
this also holds. 3 For the vector and axial tastes, however, the relationship is much more 
complicated in the non-degenerate case, although we will not need the explicit relationship 
in that case. 

Since we are assuming that the axial rf mass-squared is going negative, we are going to 
keep only this meson in $. This is valid right near the critical point and we assert that this 
remains valid so long as no other mass-squared is negative in the new phase. If this occurs, 



2 We remark that often in continuum XPT this relationship is used between the two bases. In that case, 
when m u ^ md, there are additional mass- mixing terms between the ir° and the rj which are usually 
neglected. 

3 We note that we need not write the tensor and pseudoscalar taste mesons in the physical basis, since they 
do not mix. Since this is an orthogonal transformation, the bases are equivalent. Writing them in the 
physical basis makes the subsequent discussion simpler. 
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then there will be a more complicated pattern of symmetry breaking. We will see that this 
is a possibility, and thus the phase diagram is richer than discussed here. However, we will 
only focus on one region of the phase diagram here. 

Keeping only the r)' A , $ and £ are block-diagonal in flavor space. We can write the 
condensate (£) in terms of the five real numbers A and B^: 



{T>)=A{1 



12x12 



(17) 



with the condition that A 2 + B^B^ 
get (dropping constants) 



1. We write A = cos 9 and substitute this into W to 



W 



-3/ 2 



2fim cos 9 + 



1 



a 2 A A + ^a 2 5' A ) cos 2 9 



(18) 



To find the condensate, we wish to minimize this potential. If 5' A is greater than 5' A crit = S' c , 
then the minimum always lies at 9 = if fim > or 9 = it for fim < 0. In this phase, 
the unbroken phase, (£) = sign(/im)l 12x i2 and the mass spectrum is exactly that discussed 
above (for fim < make the replacement /im — > \fim\ p|). 

The more interesting case is when we take 5' A < 5' c . Since 5' A is negative in this region, 



we write 5' A 



\S' A \. Here, we find the potential is minimized when 



cos 9 



■ijim 



3a 2 \5' A \ -Aa 2 A A 



(19) 



so the condensate has the form 

/ 

<£> = 



exp [iOn^i^)] 

exp [iOn^it^)] 










exp [1^(1^)] j 



where is a unit vector. 

The condensate breaks the SO (4) symmetry down to SO(3), giving us three Goldstone 
bosons. As noted in the introduction, these are not exactly Goldstone bosons, because we 
are not breaking an exact symmetry. The SO (4) symmetry will be broken down to the 
lattice subgroup by terms of C(a 4 ). Thus the masses of these pseudo- Goldstone Bosons are 



m GB ~ fl4 Ar 



QCD- 



Since the condensate is proportional to the identity in flavor space, the flavor-charged 
mesons (the tt + , K°, etc.) have the same form as the 7r° and rj in the broken phase. 
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Thus, below we only list the masses for the 7r° of the various tastes, but keep in mind that 
m 2 = Tn 2 + = m 2 = m 2 K+ = m 2 K0 for all tastes. Additionally, this form of the condensate 
implies that there is no flavor mixing between the 7r°, r\ and rf in this case. We pick the 
4-direction for the condensate, so that = 5^ to get explicit results for the masses. We 
also define the quantities 

16 

S' 2A = j^C 2A , (20) 
16 

44 = J 2 C,A, (21) 

so S' A = S' 2A — 5' 5A , an d similarly for A — > V. We will organize the masses according to taste, 
by whether or not they mix. 



A. Tastes which mix 



The £4 & £5 tastes mix with each other, as do the ^4 & ^5 tastes (with i = 1, 2, 3). For 
the £4 & £5 tastes, we have the same mixing for any of the mesons. With being 7r + , 7r°, 
K + , K°, 77, or r)', we define the new fields 

0f w = cos #04 -sin #0 5 (22) 
0£ ew = sin #04 + cos #0 5 . (23) 

These new fields are defined so that at cos# = 1, the critical point, 04 6w is purely 04 and 
similarly 0g ew = 05. This means that in the chiral limit, at cos# = 0, this is reversed: 02 ew 
is purely 05 and vice versa. In terms of these new mass eigenstates, we find the masses 

a 2 A v -a 2 A A + ^a 2 \5' A \ (24) 
3 

mJ,n C w = -a 2 A A + -a 2 \5' A \ (25) 

a 2 A v -a 2 A A + ^a 2 \5' A \-^a 2 5' v (26) 

Note the combination — 4a 2 A A + 3a 2 |5^|, which appears frequently, is positive, so Eqs. (j2H) 
through ()26j) are all positive in the broken phase. 

As for the £ i4 & £ i5 tastes, we have the same mixing as in Eq. (j2^jl . with 4 — > i4 and 



m I), new 
"4 



m _(J,ncw 



new 
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5 — > i5, again with <p referring to any of the six mesons. These have the masses 

a 2 A T -a 2 A A + ^a 2 \5' A \ (27) 

(28) 
(29) 

Again, all of these squared masses are positive for our range of parameter values. We see 
that the rj' i5 new are the three expected Goldstone bosons. 



O.ncw 


ff^irJ new 


fTh O.new 


= ~/Wa 


fTb^J new 
'H5 


= 0. 



B. Tastes which do not mix 

The rest of the tastes do not mix with each other, and so we just list the masses for each 
of these states. 



3 64m 2 u 2 (a 2 A 7 - a 2 A^ 

7« \ s a\ + — T 

4 Al (Aa 2 A A -3a 2 \5' A \) 2 



<o = a 2 A V - a 2 A A + °-a 2 \5' A \ + T/ „ o. " .tT (30) 



9 9* 3 . 64m 2 u 2 a 2 A A .„ . 

m 2 = - a 2 A A + -a 2 5' A + ^— * 31 

4 1 Al (4a 2 A A -3a 2 |<^|) 2 

(a 

4" 1-/11 ' (4a 2 A7-3a 2 |^| 



< s = + 7 a 2 |^| + ^ (32) 

9 9 a 9 a 3 . 64m 2 u 2 (a 2 Ay - a 2 A T ) 
m 2 = a Ay - a 2 A A + -a 2 |^| + ^ v f— r; (33) 



m 



9a 9a 3 9|P . . 16m u (4a Ay — 4a Ay + 3a 

2 , = a 2 A T - a A^ + -a 2 (ft + — =-= — 

" l 4 1 Al (4a 2 A A -3a 2 |<^|) 2 



m 2 



3 2 64m 2 /i 2 (a 2 Ay - a 2 A A ) 

4 a A ' + (4a 2 A^-3a 2 |(ftJ) 



a 2 |^| + :/ J:", (35) 



2 3 2 1 f/ I ^ Sri 2 r/ 



16m? (i? (4a 2 Ay - 4a 2 A A + 3a 2 5' 2V + 3a?6' 5V , 
+ (4a 2 A A -3a 2 |^|) 2 ' 

In the chiral limit, when m = 0, we have several degeneracies, similar to those found in 

Ref. In this limit, the states which mix are purely one state (e.g., 7T^ CW = 7if 5 ), and 

we denote these by the pure state. When m — 0, the U(1) A is restored, and we see that 

m 2 = m 2 = m 2 ,, m 2 = m 2 = m 2 , and m 2 = m 2 . Additionally, we have m 2 = m 2 , 

w iB n i "i W 4 w 45 '45 w 5 n I ' 7r i4 n ij 

which occurs in the massless case due to an accidental 5*0(3) symmetry about the minimum 
of the potential when 5' A = 0. This degeneracy recurs at the critical point for a different 
reason: due to the restoration of the 5*0(4) symmetry. Note that the limits of m — > and 
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5' A — > o^. do not commute. If we take 5' A — > 5£ first, the masses all become those of the 
unbroken phase, so taking m — > gives the massless forms of Eq. (jHJ). However, taking 
m — > first while in the broken phase then the 5' A — > 5' c limit results in a much different 
spectrum which is not the same as in the unbroken phase. 

The appearance of the "incomplete" hairpin parameters S' 2V and S' 5V in the rj'^ mass shows 
that the phase structure is more complicated than just the two phases we have discussed. 
Since we do not know the size or the sign of either of these two parameters, these could 
cause the mass-squared of the 7^- to go negative. If this is the case, the spectrum we have 
laid out will not be strictly valid, and we would have to take into account the possibility of 
a condensate also in the r/'-directions. 

We can determine the requirements on the parameters to cause the rj'- mass-squared to 
be negative. This happens when 



So in the massless case (where we will see later that we are most likely to be in the broken 
phase), this condition requires a 2 5' 5V + a 2 5' 2V > a 2 \5' A \. So the further we go into the broken 
phase, \5' A \ gets larger, meaning that we have a less chance of finding m 2 , < 0. Here we are 



studying the simplest case, where this criterion is not met, and the spectrum above is valid. 

We cannot know the values of these parameters from first principles; we must calculate 
them from fits to lattice data. This can be done by calculating quantities for non-Goldstone 
mesons (such as the taste-singlet or taste- vector charged pions). The specific combination 
5 A = 5' 2A — S' 5A (and similarly for A — > V) appears when calculating the pseudoscalar taste 
masses and decay constants, but in other taste mesons, we find that S' 2A +5' 5A also appears. In 
order to have a complete determination of the taste- violating parameters, one must calculate 
the expressions for the masses and decay constants for each of the other tastes and match 
those expressions to lattice simulation data. 

This will not change our results in Sec. however, since those results are focused on 
studying the location of the critical point, S' c . The criterion 5' A > 5' c is required for any 
symmetry breaking; if this does not hold we have the mass spectrum of the unbroken phase. 
If Eq. (pTTj) holds, this would change the spectrum of masses and the breaking pattern 
as well as the condensate itself. We will discuss the functional dependence of S' c on the 
various parameters, and examine the possibility of current (and near future) simulations 



a 2 S' 5V + a 2 S' 2V > 



3a 2 \5' A \(3a 2 \5' A \ - 4a 2 A A ) 2 + 256(a 2 A T - a 



2 A A )id 2 m 2 



(37) 



3(3a 2 |^| -4a 2 A^) 2 - 192/i 2 m 2 
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being performed in the broken phase. If we are in the broken phase, then we would have 
to determine if Eq. (|3*7)l holds; if not, the spectrum discussed above would be correct. If it 
does, a more complete study should be done of the resulting mass spectrum. As it is, a very- 
rich phase structure is possible for staggered fermions. 



IV. THE BROKEN PHASE FOR 2+1 FLAVORS 



In the 2 + 1 case, we set m u = m^ = mi ^ m s , so M. = diag(m;7, mil, m s I). Here when 
diagonalizing the flavor neutral axial taste mass matrix, the mass eigenstates for the axial 



tastes to be given by 



m 2 

A 



2\imi + a 2 A A 
1 



m 



rn 



>IA 
2 



Z = 



2/j(mi + m s ) + 2a 2 A A -a 2 5' A + Z 



2ji{mi + m s ) + 2a 2 A A + -^l 2 5 a — Z 



4/i 2 (m s - mi) 2 - a 2 5' A n (m s - mi) + 



(38) 
(39) 



The vector tastes have the same form (with V — > A), but again we are only focusing on the 
axial taste masses vanishing. 

0ur definitions of „ a nd rf A a r e different from that in Ref. Q. We define the nA a, tfie 



field which becomes degenerate with the n A when mi = m s . Thus the tja from Ref. Q| is 
now our r]' A since 5' A < (as determined from simulations). Using this definition, it is the r)' A 
mass which vanishes at the critical value of 6 A , just as in the 3-flavor case. In the 2+1-flavor 
case, this critical point is given by the expression 



a% = - 

The condensate has the form 
/ 



4 (2/im; + a 2 A A ) (2{xm s + a 2 A, 
2/x (mi + 2m s ) + 3a 2 A A 



(40) 



exp [iOn^ii^) 











exp [i0n^(i^ 5 )] 
V u exp [mOn^i^s)] J 

n M is a unit vector in the \i direction and we have defined 

8/i(m s — mi) — a 2 5 



(41) 



* 2J )' A - AZ 



2a 2 5' A 



(42) 
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which arises from the relationship between the rf A and the flavor basis 

t/ a =M(U a + D a + kS a ) , (43) 

with Af a normalization constant. Note that in the limit that 5' A — > — oo, or m s — > mi, 
k, — > 1, giving us the same definition of the f]' A as in the 3-flavor case. 

9 in Eq. (}4*Tj) is the solution which minimizes the potential, which in this case has the 
form 



w = -/ 



4/im/ cos(6>) + 2yum s cos(k0) + I a 2 A A + -a 2 5^ ) cos 2 (6 l ) 



2 



+ i ( a 2 A A + -a 2 5' A J cos 2 (k0) - -a 2 ^sin(6») sin(/«0) 
2 V 4 / 2 



+ const. (44) 



We plot this potential in Fig. [T]for the unbroken phase (so that > In this and 
subsequent figures, we are plotting the dimensionless quantity W = a 2 W/ f 2 as a function 
of 9. There are two minima, one of which is the global minimum at 9 = 0, which defines 
the vacuum state as we expect in the unbroken phase. At the critical point, 5 A = 5' c , shown 
in Fig. |21 both minima have flattened out; if we were to expand Eq. f)44j) about 9 = 0, we 
would see the 9 2 term vanish here. Once we have entered the broken phase, in Fig. El the 
minima have both shifted towards each other, and the global minimum is at a non-zero value 
of 9. As 5' A — > — oo, we find that the minima correspond to two equivalent solutions for the 
ground state here. In terms of the potential and the critical point, the 5' A — ► — oo limit is 
similar to the chiral limit [see Eq. (|40|) ]. where the U(1) A symmetry is restored. In the chiral 
limit, there are two minima given by 9 = or 7T, and they are equivalent. Shifting 9 by it 
merely sends S — > — E, which is a symmetry of the Lagrangian in the massless case. In the 
case of S A — > — oo, the minima are not at or it, but roughly near 7r/2 and 37r/2, and are 
equivalent minima in the potential. 

We cannot solve for an analytic solution for cos# in this case as we did in the 3-flavor 
case, although we can solve for the masses of the mesons in the broken phase in terms of 9, 
or more appropriately cos(#) and cos(k9). As in the 3-flavor case, there is taste mixing, but 
now the condensate is not proportional to the identity in flavor space. This implies that we 
should find mixing between the rj and rj' mesons. Since the condensate is still diagonal in 
flavor space, there will be no mixing between the charged fields. The up and down quarks 
are degenerate, so the ir°, tc + and n~ have the same general form; there are only mixings 
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among the tastes but there is no flavor mixing. The expressions are complicated and not 
very enlightening, so instead of looking directly at the mass spectrum, we will study this 
phase numerically. In the next section, we will look at the critical point as a function of the 
various parameters and see the criteria for entering this staggered Aoki phase. 



V. NUMERICAL RESULTS 



With current MILC simulations, we now discuss the dependence of the critical point on 
the various parameters in our theory. We will look at data from the "coarse" (a ~ 0.125fm) 
and "fine" (a w 0.09fm) MILC lattices 0,0]. We will look first at the quark mass dependence 
on the two different lattices with all other parameters fixed, and then see what happens at 
fixed quark mass when we decrease the size of the splitting, A^. This would correspond 
to staggered actions that have smaller discretization errors than the Asqtad action, such as 
HYP staggered fermions j^, I2]] . 

Figures Eland El show the critical S' A , Eq. (f4T)|) . on the coarse and fine lattices, respectively, 
as a function of light quark mass at fixed strange quark mass. In lattice units, the strange 
quark mass is 0.05(0.031) on the coarse (fine) lattice. As noted in Ref. these values are 
slightly larger than the physical strange quark mass, but there is no qualitative difference 
in the results here when we lower am s by w 20%. The horizontal line shows the value of 5 A 
which was determined in chiral fits to the mass and decay constant of the pseudoscalar-taste 
mesons We can see that we are far outside the broken phase here, for both the 

coarse and the fine lattices, even as we take the chiral limit. The relevant points in the 
parameter space are when ami > 0.005(0.0062), since this is the range of values simulated 
currently on the coarse (fine) MILC lattices, and ami ~ 0.002(0.001), the physical value of 
the average up and down quark masses on the coarse (fine) MILC lattices. 4 

We move closer to the broken phase when all three quark masses are light: the 3-flavor 
case in the chiral limit. In this case, 5' c = —4A A /3, which is surprisingly close to the value 
of 5' A currently determined from simulations jllj. Figures |H1 and [7| shows S' c as a function 
of the strange quark mass at mi = for the coarse and fine lattices respectively. The solid 
line shows the central value of 5' A that has been determined, while the dotted lines give the 



4 These values for the physical light quark mass are found relative to the "nominal" strange quark mass, 
am s = 0.05(0.031) on the coarse(fine) lattices, using the ratio m s Jmi ~ 27 [22] ]. 
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statistical error. We can see that there is a better chance to enter the broken phase for very 
light strange quark mass, although the error is large enough that this may not occur. A 
more precise determination of the hairpin parameter would allow for a concrete statement 
about the broken phase at light quark masses. 

As is apparent from Figs. El and as a becomes smaller, we are somewhat less likely 
to be in the broken phase. This seems to imply that as a — > 0, there is little chance we 
will be in the broken phase, even close to the chiral limit. This should be expected, as this 
phase is a lattice artifact and thus vanishes in the continuum limit (note that as a — > 0, 
5' c — > — oo). This also makes it clear that at larger lattice spacing, 5' A is more likely to be 
in this phase. We anticipated this, since the broken phase appears when the discretization 
errors are comparable to the lattice quark masses. This implies that we must not simulate 
at too coarse a lattice spacing, and the values used here are sufficiently fine so as to not 
approach the broken phase. At these lattice spacings, it is clear that there is only a problem 
if all three quarks are very near zero, which is not a physical point in the parameter space. 

One may wonder how this changes with more improved actions, such as HYP staggered 
fermions |2Qj. Let us define A A as the axial-taste splitting for some more- improved action, 
and we define the quantity 

* - ^, (45) 

with the Asqtad splitting. We plot in Fig. |H] for the coarse lattice spacing a 4 S' c as a 
function of x. We only plot the coarse lattice here, since we expect this to be the place 
where the problem would be more severe. These are plotted with the lightest light quark 
mass from the MILC simulations (mi/m s = 0.005/0.05). As x is reduced, we see for both the 
coarse and fine lattices that S' c is also reduced in magnitude. For reference, HYP fermions 
lie roughly at x ~ 0.5 j2l|. Whether or not \S' A \ is reduced by the same magnitude as A^ is 
uncertain, and this is something that must be checked in a dynamical simulation. A good 
improvement program should improve all the lattice artifacts to some degree, so we would 
expect \5' A \ to be reduced as well. As can be see from the figure, as long as \S' A \ does not 
increase by a significant factor after improvement (the central value for 5 A in the Asqtad 
case is shown by the solid line), then we will not be in the broken phase. 
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VI. CONCLUSION 



We have discussed the phase structure of a lattice theory with staggered quarks, which 
has an approximate [to 0(a 2 )] SO (4) taste symmetry, as seen in the chiral theory This taste 
symmetry can be spontaneously broken down to 5*0(3) for certain values of the parameters 
in the theory. This most likely will not occur for current staggered actions, even among 
improved theories, unless all three quark masses are close to the chiral limit. At that point, 
there could be a phase transition with 3 pseudo-Goldstone bosons [with masses ~ (9(a 4 )], 
and a rather different mass spectrum than in the unbroken phase. 

The possibility of a flavor-symmetry (or in this case taste-symmetry) breaking phase is 
a general feature of any lattice theory which has additional terms in the squared meson 
masses at finite lattice spacing. While originally discussed in the context of Wilson fermions 
(Aoki phase), we see that there are complicated ways for the taste symmetry with staggered 
fermions to be broken as well. The difference here is mainly that the taste symmetry is 
already broken by the lattice spacing, and it is only the approximate 5*0(4) symmetry 
which is broken. So unlike the Wilson case, there are significant corrections at higher order 
in the lattice spacing. 

This brings up an interesting question about the 0(a 4 ) corrections. We are stating here 
that when the lattice spacing becomes too large, we are most likely to enter this broken 
phase. However, as a increases, higher-order scaling violations become more important. It 
is unclear whether or not the broken phase will remain at extremely coarse lattice spacing. 
Using a chiral Lagrangian that includes higher order corrections |23j], this could be studied 
by determining the relative sizes of higher-order contributions. 

An interesting question arises relating to the quenched approximation of staggered 
fermions. For Wilson fermions, there is the possibility of entering the Aoki phase, and 
this can cause significant problems with locality when simulating Wilson, Domain Wall, or 
Overlap quarks 24]. For the staggered case, however, there is no Aoki phase in this context. 
There cannot be an Aoki phase in the context of the current analysis, since the S' A splittings 
do not add to the mass of the rj' A . So unlike the case of Wilson quarks, there is no difficulty 
when simulating staggered quarks in the quenched approximation. 
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FIG. 1: Potential in the unbroken phase; parameters here are for the coarse MILC lattice, with 
the light/strange lattice quark masses 0.005/0.05. The global minimum is clearly at = 0. 
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FIG. 3: The potential in the broken phase, where the global minimum now is at a non-zero value 
of 9. Note the other (local) minimum has shifted down and to the left; as S' A — > — oo, the two 
minima will both be global minima at different values of 0. 
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FIG. 4: The dashed line shows S' c as a function of the light quark mass (in lattice units) on the 
coarse lattice, for a lattice strange quark mass of 0.05. Th e solid horizontal line is the value of 5' A , 



determined from chiral fits to the 2 + 1 MILC data 
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FIG. 5: Same as Fig. |1] but for the fine lattice spacing, and for a lattice strange quark mass of 
0.031. 
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FIG. 6: The dashed line is S' c on the coarse lattice as a function of the strange quark mass for 
mi = 0. The solid line is the value of 5' A , with the statistical errors given by the dotted lines. For 
very light m s we see that there is a significant probability that we will enter the broken phase. 
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FIG. 7: Same as Fig. El but for the fine lattice spacing. There is still a significant probability that 
we could be in the broken phase for light m s but it is somewhat less likely than on the coarse latice. 
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FIG. 8: S' c versus the fractional improvement of the splitting for the coarse lattice. _x = 1 
corresponds to the Asqtad action, while x ~ 0.5 corresponds to HYP staggered fermions 
The horizontal solid line is the value of 5' A for the Asqtad action for reference. 
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